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Abstract: We consider one- loop amplitudes of a Higgs boson coupled to gluons in the 
limit of a large top quark mass. We treat the Higgs as the real part of a complex field cj> 
that couples to the self-dual field strengths and compute the one-loop corrections to the (fi- 
NMHV amplitude, which contains one gluon of positive helicity whilst the remaining three 
have negative helicity. We use four-dimensional unitarity to construct the cut-containing 
contributions and a hybrid of Feynman diagram and recursive based techniques to deter- 
mine the rational piece. Knowledge of the 0-NMHV contribution completes the analytic 
calculation of the Higgs plus four gluon amplitude. For completeness we also include ex- 
pressions for the remaining helicity configurations which have been calculated elsewhere. 
These amplitudes are relevant for Higgs plus jet production via gluon fusion in the limit 
where the top quark is large compared to all other scales in the problem. 
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1. Introduction 



The search for the Higgs boson is one of the primary objectives of the Large Hadron Collider 
(LHC). If discovered, attention will swiftly turn to the exploration of the Higgs sector and 
the measurement of the Higgs couplings to the weak gauge bosons and to fermions. The 
main Higgs production processes at the LHC are gluon fusion, which proceeds via a top 
quark loop, and vector boson fusion (VBF), which is dominated by the t-channel exchange 
of weak bosons. The Higgs production rate via vector boson fusion is typically about an 
order of magnitude smaller than the gluon fusion rate, but has a characteristic signature 
with two forward quark jets. The next-to-leading order (NLO) QCD corrections to the 
VBF process are small (5-10%) [1-4] which makes vector boson fusion very attractive for 
the measurement of the Higgs coupling to weak bosons. Recently, the full EW+QCD 
corrections to this process have been computed [5, 6]. 

The dominant background to VBF comes from Higgs plus two jet production via gluon 
fusion. Leading order calculations have been performed both with the exact mt dependence 
[7] and in the large mt limit [8, 9] where the top quark loop is replaced by an effective local 
interaction C (mt)HG^ G ^ u [10-12]. For inclusive Higgs production it has been shown that 
this approximation is valid over a large range of Higgs masses [13]. The approximation 
remains valid for processes with increased numbers of jets provided that the transverse 
momenta of the jets is smaller than m t [14]. 

The next-to-leading order QCD corrections to the gluon fusion Higgs plus two jet rate 
have been calculated numerically in the large-m^ limit [15] using the real radiative correc- 
tions of Ref. [16-18] and the partly analytic, partly numerical virtual one-loop corrections 
of Ref. [19]. There are four distinct processes that contribute to the virtual NLO correc- 
tions, gg — ► ggH, gg — ► qqH, qq — ► qqH , qq — > q'q'H and the associated crossings. Ref. [19] 
provided analytic expressions for the spin- and colour-averaged one-loop "squared" matrix 
elements for qq — > qqH and qq — > q'q'H and evaluated the other two processes numerically. 
Because of the size of the NLO QCD corrections, for studies of Higgs phenomenology it is 
important to incorporate the NLO QCD corrections to Higgs plus two jets in an efficient 
and flexible way. This has led to a collective effort to derive compact analytic expressions 
for the one-loop corrections for the gg — > ggH and gg — ► qqH processes. 

An effective way of deriving compact analytic expressions for the Higgs plus four parton 
amplitudes is to employ on-shell unitarity methods. The original unitarity method involved 
using four-dimensional double cuts of amplitudes to classify the coefficients of discontinu- 
ities associated with physical invariants [20, 21]. Unitarity cuts found many applications in 
the calculations of amplitudes in N = 4 and M = 1 supersymmetric Yang-Mills theories, 
since these amplitudes can be fully constructed from their unitarity cuts. Amplitudes in 
non-supersymmetric theories contain rational terms which cannot be determined from four 
dimensional cuts but can be computed using a number of complementary methods [22-27]. 

In recent years the unitarity method has been generalised to include multiple cuts [28- 
35] , this has lead to an explosion in the number of unitarity based techniques in calculations 
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of NLO processes. In particular by working in D-dimensions [22, 36-39], one can completely 
determine the amplitude from its unitarity cuts. This is because the rational parts of 
one-loop amplitudes arise from higher-dimensional basis integrals, which are sensitive to 
L>-dimensional cuts. At tree-level the discovery of another on-shell method, the BCFW 
recursion relations [40, 41], also sparked developments at one-loop. Recognising that the 
rational part of one-loop amplitudes also satisfied recursion relations allowed a fully four 
dimensional on-shell method to be introduced (the unitarity bootstrap) [23, 42]. Several 
groups have produced sophisticated numerical programs based on on-shell methods which 
aim to efficiently calculate NLO one- loop amplitudes of relevance to the LHC [43-45]. 

In order to efficiently compute one-loop amplitudes using on-shell unitarity based tech- 
niques, it is desirable to have compact tree amplitudes. For amplitudes involving a Higgs 
boson coupling to partons, it is convenient to split the real Higgs scalar into two com- 
plex scalars ((ft and (ft) such that H = (ft + (ft [18]. In this case, the effective Higgs- 
gluon interaction also separates into two parts - the (ft couples directly to the self-dual 
gluon field strengths, whilst (ft couples to the anti-self-dual gluon field strengths. The 
tree ^-amplitudes are MHV-like and have a very simple structure. <ft amplitudes are ob- 
tained from (ft amplitudes by complex conjugation. With this breakdown, amplitudes for 
a light pseudo-scalar Higgs boson A can be constructed from the difference of (ft- and <ft- 
amplitudes, A ~ (ft — (ft . The ^4-amplitudes may be relevant amplitudes for SUSY theories, 
provided that tua < 2mj. 

Analytic calculations of one-loop (^-amplitudes to date have not included NMHV he- 
licity configurations. Previous one-loop calculations for amplitudes with (ft coupling to 
any number of gluons have restricted the number of negative helicity gluons to: none or 
one (which are rational amplitudes) [46], two (the (ft-MHV configurations) [47, 48]. An 
expression for the amplitude containing only negative helicity gluons has also been ob- 
tained [49]. Recently new analytic results for the helicity amplitudes of the (ftqqg ± g T and 
(ftqqQQ processes have appeared [50]. 

In this paper, we calculate the one- loop NMHV amplitude A^\(ft, 1+, 2~ , 3~, 4~) and 
combine it with the already known A^\<ft , l + ,2~,3~,4~) amplitude [46] to complete the 
set of analytic helicity amplitudes for the process — > Hgggg. The only remaining missing 
analytic piece for the process pp —> H + 2j is the (ftqqgg-NMWV amplitude. 1 

By employing a four-dimensional unitarity-based strategy, we reconstruct the cut- 
constructible term as a combination of n-point scalar integrals (with n = 2, 3, 4) which are 
referred to as bubble-, triangle- and box-functions. By applying a specific set of cuts both 
to the amplitude and to the basis integrals, one can isolate specific coefficients appearing 
in the basis expansion. Cutting a propagator means restricting the momenta so that 
the propagating particle is on-shell, and for four dimensional loop momenta, setting n 
propagators on-shell leaves 4 — n of the loop momentum free. Hence, in the case of the 
quadruple-cuts [28], the loop momentum is completely fixed, as result, the determination 
of the box-coefficients is reduced to an algebraic operation. To compute triangle- and 

x We note that the <pqqgg-NMHV amplitude has recently been computed [51]. 
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bubble-coefficients, fewer cuts must be applied. In these cases the free components of the 
loop momentum become integration variables of phase-space integrals which can be carried 
out with the mathematical methods of complex analysis. In our calculation we determine 
the coefficients of the triangle- and bubble-functions using two variants of the triple- and 
double-cut integration techniques [29, 30, 52]. 

The rational part of the amplitude however, cannot be detected by four-dimensional 
cuts and as a result some other method must be used. We find it convenient to separate 
the rational part into two terms, one being sensitive to the number of active light flavours 
and one which is not. The former piece is efficiently computed using Feynman diagrams, 
whilst the latter can be derived using on-shell recurrence relations [40, 41]. 

This paper is organised as follows, in section || we describe the model for Higgs inter- 
actions in the large m% limit. In particular we describe the breakdown of Higgs amplitudes 
into cj) and $ contributions and discuss the colour decomposition of one-loop Higgs plus 
multi-gluon amplitudes. We separate the one-loop primitive amplitude into two parts: C4 
which contains the cut-constructible parts of the amplitude, and R4 which contains the 
remaining rational pieces. In section |3] we use various unitarity methods to determine the 
coefficients of the various one-loop basis integrals which appear in the cut-constructible 
C4 and we derive the rational term R4 in section |j. In section [5| we present all the prim- 
itive helicity amplitudes for the — > Hgggg process. Section |6| contains numerical values 
for the Higgs helicity amplitudes at a particular phase space point. We draw our conclu- 
sions in section |7[ We also enclose two appendices containing the tree expressions used as 
inputs into the unitarity based calculations and the definitions of the finite parts of the 
one- loop basis integrals. Throughout our calculations we have made extensive use of the 
S@M package [53]. 



2. The Model For Higgs Interactions In The Large m t Limit 



We work in an effective theory where the Higgs couples to gluons through a top quark loop, 
but the top quark degrees of freedom have been integrated out, 

£ mt = C HtTG ^Q» r (2.1) 

The effective coupling C has been calculated up to order C(af) in [54]. However, for our 
purposes we need it only up to order 0(a 2 s ) [55], 

c = ^( 1 + ^) +0(a ' ) (2 ' 2) 

where v = 246 GeV and the strong coupling constant is a s = g 2 /(4?r). 

Following Ref. [18] we introduce the complex field (j) (and its conjugate (jy 

, (H + iA) ,, (H-iA) , , 

0= \ L , f = ; - (2-3) 
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The effective interaction linking gluons and scalar fields also splits into a piece containing 
<f> and the self-dual gluon field strengths and another part linking (fy to the anti-self-dual 
gluon field strengths. The last step conveniently embeds the Higgs interaction within the 
MHV structure of the QCD amplitudes. The full Higgs amplitudes are then written as a 
sum of the <fi (self-dual) and <jy (anti-self-dual) components, 

4° (#;{?>*}) = 4°(6{p*}) +4 ) (0 t ,{»})- (2.4) 

We can also generate pseudo-scalar amplitudes from the difference of 4> and (jr components, 

4 (4{p*» = 7 (4°(<M^}) -4V,fe})) • ( 2 -5) 

Furthermore parity relates (f> and <jr amplitudes, 

4 m ) (tf,g$\...,gt)= (4 m) 5r Al , • • • , 9n K )) * • (2-6) 

From now on, we will only consider ^-amplitudes, knowing that all others can be obtained 
using eqs. (|g)-(gj). 

The tree level amplitudes linking a <j> with n gluons can be decomposed into colour 
ordered amplitudes as [8, 16], 

A^{^{k l ,X h a i })=iCg' 1 - 2 £ tr(T<MD . . . J^M) 4°>(6 <x(l Al , ..,n A «)). (2.7) 

Here S n /Z n is the group of non-cyclic permutations on n symbols, and j Xj labels the 
momentum pj and helicity Xj of the j th gluon, which carries the adjoint representation 
index a*. The T ai are fundamental representation SU(iV c ) colour matrices, normalised so 
that Tr (T a T b ) = 5 ab . Tree- level amplitudes with a single quark-antiquark pair can be 
decomposed into colour-ordered amplitudes as follows, 

4 O) (0, {JH, Ai, fli}, {pj,Xj,ij}) (2.8) 
= iCg n ~ 2 £ (™ 2) • • •T a -^)i lin M4>, lV(2 A2 , • • • , in - l) A "- 1 ),n" A ) , 

<TG5 n -2 

where S n -2 is the set of permutations of (n — 2) gluons. Quarks are characterised with 
fundamental colour label ij and helicity Xj for j = 1, n. By current conservation, the quark 
and antiquark helicities are related such that Ai = — A n = A where A = 

The one-loop amplitudes which are the main subject of this paper follow the same 
colour ordering as the pure QCD amplitudes [20] and can be decomposed as [46-49], 

[n/2]+l 

A^(^{k l ,Ka t })=iCg n £ £ G n , c {a)A^ {^a{l x \ . . . ,n A ")) (2.9) 

where 

G n; i(l) = iV c tr(T ai ...T a ") (2.10) 
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G n;c (l) = tr(T ai ■■■r 0c - 1 )tr(T 0c •••T a ") , c> 2. (2.11) 

The sub-leading terms can be computed by summing over various permutations of the 
leading colour amplitudes [20]. 

In this paper we define kinematic invariants associated with sums of gluon momenta 
as follows, 

Sij = {Pi +Pj) 2 , s ijk = (pi + pj + pk) 2 etc. (2.12) 

We will express helicity amplitudes using the notation of the spinor-helicity formalism, 

(ij) =u_(k i )u+(k j ), (2.13) 
[ij] =u+(k i )u-(k j ), (2.14) 

where u±(ki) represents a massless Dirac spinor associated with either positive or negative 
helicity (and a momentum fcj). Spinor products are related to kinematic invariants through 
the following relation, 

Sij = (2.15) 
Chains of spinor products are written as 

(i\j\k] = (ij)[jk] (i\jk\l) = mjk}(M), etc. (2.16) 
For example, using momentum conservation we have, 

4 

<iMfc]=-£fom (2-17) 
i=i 

Throughout this paper we use the following expression for the </>-NMHV tree amplitude 
4°)(<M + ,2-,3-,4-) = 

m iW + (jMil! (2Mi] 3 (218) 



si24<12)(14)<2|p^|3]<4|p^|3] si23<4|p^|3][12][23] s 134 (2|^|3] [14] [34] ' 

This compact form can be derived using the BCFW recursion relations [40, 41] and agrees 
numerically with the previously known expression derived from MHV rules [18]. It clearly 
possess the correct symmetry properties under the exchange {2 <-> 4}, and factors onto the 
correct gluon tree amplitude (which is zero) in the limit of vanishing p^. The other tree 
amplitudes we require for this work are listed in Appendix A. 



3. Cut-Constructible Contributions 

We choose to expand the one-loop primitive amplitude in the following form, 

4% 1+, 2", 3-,4") = c r (C 4 (0, 1+ 2-, 3", 4") + R A (<t>, 1 + , 2", 3" 4")), (3.1) 



- 5 - 



where 

r 2 (i- e )r(i + 6) 

Cr = (An)^T(l - 2c) ■ (3 ' 2) 

In fl3.l|) , C4(0, 1 + , 2~, 3~, 4~) denotes the cut-constructible parts of the amplitude, whilst 
Ri(<f>, l + ,2~,3~,4~) contains the remaining rational pieces. In this section, we focus our 
attention on C^cfr, 1 + , 2~, 3~, 4~), while an analytic expression for R±{4>, 1 + , 2~ , 3~, 4~) is 
derived in section ||. 

We employ the generalised unitarity method [28-32, 52] to calculate the cut-constructible 
parts of the one-loop amplitude. We can further decompose C4 in ( p.l| ) into a sum over 
constituent basis integrals, 

C7 4 (</»,l+,2-,3-,4-) =^C7 4;J J 4;i + ^C3 ;4 X3;i + ^C 2;i J 2;4 . (3.3) 

iii 

Here Ij-i represents a j-point scalar basis integral, with a coefficient Cj ; j. The sum over i 
represents the sum over the partitions of the external momenta over the j legs of the basis 
integral. 

Multiple cuts isolate different integral functions and allow the construction of a lin- 
ear system of equations from which the coefficients can be extracted. When considering 
quadruple cuts of one-loop amplitudes, one is forced to consider complex momenta in order 
to fulfill the on-shell constraints [28]. The four on-shell constraints are sufficient to isolate 
each four-point (box) configuration by freezing the loop momentum, thereby allowing the 
determination of the corresponding coefficient by a purely algebraic operation. To isolate 
the coefficients of lower-point integrals, one needs to cut fewer than four lines. In this 
case the loop momenta is no longer completely determined, but, according to the number 
of cuts, some of its components are free variables. In this case the computation of the 
three- (triangle) and two-point (bubble) coefficients can also be reduced to algebraic pro- 
cedures by exploiting the singularity structure of amplitudes in the complex-plane, either 
by explicit subtraction [45, 56] or by Laurent expansion [29]. Alternatively one can extract 
the coefficients of bubble- and triangle-functions by employing the spinor-integration tech- 
nique, which can be applied to both double- [31, 32] and triple-cuts [30]. This method has 
recently inspired a novel technique for evaluating the double-cut phase-space integrals via 
Stokes' Theorem applied to functions of two complex-conjugated variables [52]. 



3.1 Box Integral Coefficients 

We begin our calculation of the </>-NMHV amplitude by computing the coefficients of the 
scalar boxes using generalised unitarity with complex momenta [28]. In general there are 
sixteen box topologies, which can be obtained from cyclic permutations of those shown in 
Fig. [I]. When a box graph contains a sequence of three point vertices, a non-vanishing 
solution is only found when the vertices alternate between MHV and MHV-types [28]. For 
the specific helicity configuration we consider this is not possible for the graphs of Fig [!](&); 
the three-point MHV vertex involving (f> and two loop-gluons cannot have both adjacent 
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Figure 1: The various box integral topologies that appear for j4 4 (0,1,2,3, 4). From the four 
topologies we must also include cyclic permutations of the four gluons. Here (a) has one off-shell 
leg (one-mass) whilst (b)-(d) have two off-shell legs. In (6) the two off-shell legs are not adjacent 
and we refer to this configuration to as the two-mass easy box, while in (c) and (d) the two off-shell 
legs are adjacent and we label them as two-mass hard boxes. 
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three-point gluon vertices to be of MHV-type. Therefore, the coefficients of the two-mass 
easy boxes (Fig |](o)) are an zero. 

Of the remaining 12 coefficients, a further 5 are related to each other by the {2 <-> 4} 
symmetry, 

(3.4) 
(3.5) 
(3.6) 
(3.7) 
(3.8) 



3 

Q ; ,i|2|3|4(0,l + ,2-,3-,4-) = 2{1 ^^$ mM , (3-9) 

r U 1+9- 1-4-^ ^34^4i(2|p |l] 3 ^34 g 4i(34) 3 m^ 

C 4; ,2|3|4|i(0,l ,2 ,3 ,4 ) - 2si 3 4(2 |^|3][34][41] + 2 Sl34 (l|^|2](3K|2](41)- (3 - 10) 

We also find that three of the two- mass hard boxes (Fig. ^(d)) have coefficients related to 
the coefficients of eqs. Q3.4|), ( jlO| ) and ( |3.10| ), 

Q ; <A|i2|3|4(<A,l + ,2-,3-,4-) = £^ 3 £MQ^ 4|1| 2| 3 (0,l+,2-,3-,4-), (3.11) 

' ' ' S23S12 

^01231411(0, 1 + , 2", 3", 4") = £^4£41 ^,2,314(0, 1+ , 2" , 3" , 4" ) , (3.12) 

S23S34 

C 4 ;0|34|l| 2 (0, 1 + , 2", 3", 4~) = £H4£12 C 4; ^|3|4|l(0, 1 + , 2~ , 3~ , 4~ ) . (3.13) 

S 4 1S 34 

The final two-mass hard box coefficient is, 
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Q^|3|4|i2(0,l + ,2 ,3 ,4 ) 2^(3|^|2K3|p |4][21][41] + (12)(14)(2^|3](4|^|3]J 

(3.14) 

The remaining one-mass box configuration C4 ; 03|4|i|2 is the- only one which receives contri- 
butions from Nf fermionic and N s scalar loops, 



C4;</.3|4|l|2(<M + ,2 ,3 ,4 ) = S41S12 



•S124S34 



-c 



4;<£|3|4|m 



1+2- 3- 4~) 



_ N£\ 2(3|^|1] 2 _ / _ Nf_ NA [12][41](3^|2](3| P0 |4] \ 

47VJ S124 [24P V JV e JVcj ^124[24]4 ^ [ ' ' 



Each of the coefficients ( |3.10| ), (|3.9|), ( |3.14| ) and ( p. 15 ) correctly tends to zero in the soft 
Higgs limit (p^ — > 0). 



3.2 Triangle Integral Coefficients 
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Figure 2: The various triangle integral topologies that appear for a[ 1} (</>, 1, 2, 3, 4). From the six 
topologies we must also include cyclic permutations of the four gluons. (a) has one off-shell leg, 
(b)-(e) have two off-shell legs while in (/) all legs are off-shell. 

Altogether, there are twenty-four triangle topologies, which can be obtained from cyclic 
permutations of those shown in Fig. |2[ The different topologies can be characterised by 
the number of off-shell legs. Fig. ||(a) has one off-shell leg, Figs. |2|(6)-(e) have two off-shell 
legs while for Fig. all legs are off-shell. We refer to the triangle integrals with one- and 
two-off-shell legs as one-mass and two-mass respectively. They have the following form, 
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(3.16) 
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and therefore only contribute pole pieces in e to the overall amplitude. In fact, the sole 
role of these functions is to ensure the correct infrared behaviour by combining with the 
box pieces to generate the following pole structure, 

J 4( 1 )(0,l + ,2",3-,4-) = - J 4( o )(^,l+,2-,3-,4-)^^f^^y + O(6°). (3.17) 

This relation holds for arbitrary external helicities [47-49] . We computed the coefficients of 
all one- and two-mass triangles and explicitly verified eq. ( [3.17 ). The non-trivial relation- 
ship between the triangle and box coefficients provides a strong check of our calculation. 
However, we find it more compact to present the final answer in a basis free of one- and 
two-mass triangles. That is, we choose to expand the box integral functions into divergent 
and finite pieces, combining the divergent pieces with the one- and two- mass triangles to 
form ( 3.17p and giving explicit results for the finite pieces of the box functions. 

A new feature in the 0-NMHV amplitudes is the presence of three-mass triangles, 
shown in Fig. ||(/). In previous calculations [46-50] the external gluon helicities prevented 
these contributions from occurring. 

There are four three-mass triangles, which satisfy, 

C 3;0 |34|i2(0, 1 + , 2", 3", 4") = C 3 ^| 12 |3 4 (0 5 1+ 2-, 3-, 4-) (3.18) 
C 3 - m \2 3 (^ 1 + , 2- 3" 4~) = 0^123141(0, 1 + , 2" 3-, 4-). (3.19) 

The symmetry under the exchange of gluons with momenta p 2 and p^ relates the remaining 
two coefficients, 

C 3 ;0|23|4i(0, 1 + , 2", 3", 4-) = C 3 ^| 12 | 34 (0 J 1+, 4-, 3-, 2~). (3.20) 

To compute C^^^i we use both Forde's method [29] and the spinor integration technique 
[30]. For a given triangle coefficient C s .x 1 \k 2 \k 3 { ( I ) i l + > 2~, 3 _ , 4~) with off-shell momenta 
K\, K 2 and K3, we introduce the following massless projection vectors 

1 7 i*-K\Kl ' 
" 7 J*-KlKl ' 



j±{K u K 2 ) = K 1 -K 2 ± ■ K\ - K\K\. (3.21) 
In terms of these quantities we find, 

„ m^(^2) 3 (34) 3 
CWli2|34(0, 1 + , 2", 3", 4") = V o , , : , (3-22) 

which, as expected, correctly vanishes in the soft Higgs limit {ps — > 0). 
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C2;<t> 4 4 C2 ;< /,4 ' 3 C2 ; 034 
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Figure 3: The three bubble integral topologies that appear for A± (cf>, 1, 2, 3, 4). We must also 
include cyclic permutations of the four gluons. 

3.3 Bubble Integral Coefficients 

The non-vanishing bubble topologies for the </>-NMHV amplitude are shown in Fig. |3|. 
We find that the double-cuts associated with Fig. H(a) contain only contributions from 
boxes and triangles, and therefore the coefficient of log(si234) is zero. In a similar fashion, 
the double cuts associated with diagram Fig. g(c) with two external gluons with negative 
helicity emitted from the right hand vertex have only box and triangle contributions, so 
that the coefficients of log(s23) and log(s34) are also zero. 

The leading singularity of the bubble integral is C(l/e), 

However for the total amplitude there is no overall e pole, and this implies a relation 
amongst the bubble coefficients such that, 

4 

(C2;# + C^kk+l) = 0- (3.24) 

k=\ 

It is therefore most natural to work with log's of ratios of kinematic scales (rather than 
log(s//x 2 )), since the coefficients of individual logarithms must cancel pairwise. To this 
end, we express our result in terms of the following functions, 

< 3 - 25 > 

Using the Stokes' theorem method [52], we generated compact analytic expressions for 
the coefficients of each bubble-function, which we also checked numerically with Forde's 
method [29] . The combination of all double-cuts is given by, 

l+,2-,3-,4-) _ (4- § ) C*» + (l - ^ + I) Cf> (3.26) 

with 

c« _ - { "Ml, (5121 , , 12) - (5123 , Sl2) ) _ ( (2 „ 4) 1 (3 . 27 ) 

I si 24 [42J si23[32J J [ J 
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and 
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+ V <24)[42]3 3 Sl24 [42] ) M*m,'*) 

, (3| P0 |1] (4g 124 (34)[41] + (3|p^|l](2s 14 + s 24 )) r , , 
+ Sl24 (24)[42]3 Lo (S124 ' Sl2) 

2s 123 (23)(34) 2 [31] 2 r (23)(34)[31](4|p |l] r 

^T^j L 3 I.S123, S12J H ^j^j -^2 (Sl23, «i 2 J 

In the above formulae (and those following) we stress that the symmetrising action applies 
to the entire formula, and also acts on the kinematic invariants of the basis functions. We 
see that C 2 (0, l + ,2~,3~,4~) vanishes in the soft Higgs limit — ► 0. 

3.4 The Cut-Completion terms 

The basis functions L 3 (s,t) and L 2 (s,i) are singular as s — ► t. Since this is an unphysical 
limit one expects to find some cut-predictable rational pieces which ensure the correct 
behaviour of the amplitude as these quantities approach each other. These rational pieces 
are called the cut-completion terms and are obtained by making the following replacements 



in (]3T28|) 



L 3 (s 


t)~ 


->L 3 (s,t) 


= L 3 (s 


t) 


L 2 (s 


t)~ 


^L 2 {s,t) 


= L 2 (s 


t) 


Lx(s 


*)- 




= Lx(s 


t) 


L (s 


*)- 


■* L (s,t) 


= L (s 


t) 



1 fl 1 

+ 



2{ s -t) 2 \s ' ty 
1 + - 



2(s-t)\s ty 



(3.29) 



4. Rational Terms 



We now turn our attention to the calculation of the remaining rational part of the ampli- 
tude. In general the cut-unpredictable rational part of plus gluon amplitudes contains 
two types of pieces, a homogeneous piece, which is insensitive to the number of active 
flavours and a piece proportional to (1 — Nf/N c + N s /N c ), 

R 4 (<t>, l+,2-,3-,4-) = ^(<A,l + ,2-,3-,4-) + (l-^ + ^<^(0,l + ,2-,3-,4-). 

(4.1) 
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The homogeneous term R^(<f>, 1 + , 2 ,3 ,4 ) can be simply calculated using the BCFW 
recursion relations [40, 41], 

RtW, 1+, 2-, 3-, 4-) = 2A(°\4>, 1+, 2-, 3-, 4-). (4.2) 

This contribution cancels against a similar homogeneous term for the $ amplitude when 
combining the (f> and $ amplitudes to form the Higgs amplitude. 

The Np piece allows the propagation of quarks (or scalars) in the loop, and can be 
completely reconstructed by considering only the fermion (scalar) loop contribution. Fur- 
thermore, one can extract the 4> contribution to R^ p by considering the full Higgs amplitude 
and removing the fully rational (ft contribution calculated in [46]. Since there is no direct 
Hqq coupling in the effective theory, the most complicated structure is a second-rank ten- 
sor box configuration. Of the 739 diagrams contributing to the Hgggg amplitude 2 , only 
136 contain fermion loops and are straightforward to evaluate. 

After subtracting the cut-completion and homogeneous rational terms from the explicit 
Feynman diagram calculation the following rational pieces remain. 

R N P(R 1+ 9 - fl/ (23)(34)(4|p g |l][31] (3|p H |l] 2 <24><34><3|p g |l][41] 

\2\ 3 Sl23 (12)[21][32] s 124 [42] 2 + 3 Sl24 s 12 [42] 
[12] 2 (23) 2 (24)(s 23 s 24 + s 23 s 34 + s 2i s u ) 
si 4 [42] 2 3(12) (14) [23] [34] [42] 

<2|p g |l](4| Pg |l] 2[12](23)[31] 2 \) f } 
+ 3. 234 [23][34] 3[23] 2 [41][34]J/M ( T 

The last line in the above equation is the one-loop rational expression for the (fy contribu- 
tion [46]. We can thus define the rational terms for the (j) contribution. 

l/<23)<34)<4|p H |l][31] <3|^|1] 2 , (24)<34)<3|p H |l][41] 



Ri (<M ' ' 3 ' j " \2\ 3 Sl23 (12)[21][32] + 



3si 24 s 42 [42] 

[12] 2 (23) 2 (24)(s 23 s 24 + s 23 s u + s 24 s 34 
si 4 [42] 2 ' ~~ 3(12) (14) [23] [34] [42] 



+ { (2 ~ 4) 

(4.4) 



5. Higgs plus four gluon amplitudes 



In this section we present complete expressions for the one-loop amplitudes needed to 
calculate the process — ► Hgggg at NLO. 



The one-loop amplitudes presented here are computed in the four-dimensional helicity 
scheme and are not renormalised. To perform an MS renormalisation, one should subtract 
an MS counterterm from A^~\ 




2 Feynman diagrams were generated with the aid of QGRAF [57]. 
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The Wilson coefficient ( [2.2j ) produces an additional finite contribution, 

A^^A^ + crllAf. (5.2) 

We choose to split the un-renormalised amplitude into (completed) cut-constructible 
pieces and rational terms. We also separate the infra-red divergent and finite parts of 
the amplitude. The basis functions for the finite part of the cut-constructable pieces are 
one-mass and two-mass boxes, three-mass triangles, and completed functions Lj(s,i) of 
eq. ( 3.29 ). We define the finite pieces of the box and three-mass triangle integrals in 
Appendix B. 

We express a generic helicity configuration in the following form 
A { 1 ] (H, 1 Al , 2 Aa , 3 Aa , 4 A4 ) = c r (C 4 (H, 1 Al , 2 Aa , 3 A3 , 4 A4 ) + R 4 (H, 1 Al , 2 A2 , 3 A3 , 4 A4 ) ) , 

(5.3) 

where C 4 represents the cut-constructible part of the amplitude and i? 4 the rational pieces. 
We further separate C 4 (-£f, l Al , 2 A2 , 3 Aa , 4 A4 ) into divergent and finite pieces, 

C 4 (#,l Al ,2 A2 ,3 A3 ,4 A4 ) = Va(H, l Al ,2 A2 ,3 A3 ,4 A4 ) + F 4 (#, l Al , 2 A2 , 3 As , 4 A4 ). (5.4) 

The divergent part V4 contain the e singularities generated by the box and triangle contri- 
butions, and which satisfy the helicity independent infrared singularity condition, 



1 / 4 

y 4 (#,l Al ,2 A2 ,3 A3 ,4 A4 ) = - A(°)(F,l Al ,2 A2 ,3 A3 ,4 A4 )-2 I 



/' 2 



(5.5) 



The remaining cut-constructible and rational terms are finite, and depend non-trivially on 
the helicity configuration of the gluons. 

5.1 The all-minus amplitude A^(H, 1~ , 2~ , 3~, 4") 

The all-minus amplitude is symmetric under cyclic permutations of the four gluons. The 
finite part (of the cut-constructible piece) is [49] , 

Wl-,2-,3" 4") = { - 2[12][2 ^34 ][41] (\Flr(s 123 ,s 2U ; m %,s 23 ) 

+ - Flf e (si23 ; S124! TTl%, Su) + F\^(s 2 3, S34; S 234 ] 



+ |(1 «- 4), (2 <- 3)j + |(1 <- 2), (3 «- 4) j. + |(1 <- 3), (2 <- 4) 

(5.6) 

while the rational part is given by [46, 49] 

^(i/,1-,2-,3-,4-) . (i(l - £ + %)( - + S + 2 < 34 >< 41 > 



3\ N c JVJV s 123 [12] 2 [23] 2 [12] 2 [12][23] 

+ cyclic permutations. (5-7) 



S12S34 + S123S234 - s\ 2 



2[12][23][34][41] 
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5.2 The MHV amplitude A%\h, 1~ , 2~, 3 + , 4+) 



For the MHV amplitude with adjacent negative helicity gluons there is an overall ((1 
2), (3 — ► 4)) symmetry. The finite cut-constructible part is [47], 



F 4 (i7,l-,2-,3 + ,4+) 



<12>* 



2(23)<34)(41) 



F1f C (si23,S234;"2^,S23) 



+ \ F4F C (s 23 4, S134! m 2 H , S34) + - F|p 1C (si24, S123; mjj, S12) 

+ F^p (s 23 , S34; S234) + Fjf (Si4, si 2 ; Si 24 )^ 

/ iV,\ (12) 2 [43]~ 

-H 1_ iiJ^34T Ll(Sl34 ' Sl4) 

AT/ JV fl ^ [43](13p H 2)((13p g 2) + (1432)) y 

1 " ^ + wj { 3( 134 ' 14) 



(12) 2 [43] 



^l(si34, s 14. 



3(34) 
+ {(l<-2),(3<-4) 



+ 



(1-3), (2^4) 



(5.8) 



The rational terms i? 4 have the same symmetries [47] 
# 4 (#,l-,2-,3 + ,4+) = 



iVj A^A [34] / (23)(l|p g |3] 2 <14)(3|P 12 |3] 
Nc NJ 3(34) V (34) [43] [32]s 23 4 (34) [12] [32] 



(12) 2 (12) (12)(2|P 13 |4] (12) 2 
2(34) [43] 2[12] 2[41]s 34 i 2s 4i 



+ 



(1-3), (2^4) 



+ (1-2), (3 <-4) 



(5.9) 



5.3 The MHV amplitude A%\h, 1", 2+, 3", 4+) 

The alternating helicity MHV configuration has the larger set of symmetries, (1 <-> 3), 
(2 <-> 4) and ((1 <-> 3), (2 <-> 4)). The finite cut-constructible contribution is [48], 



F 4 (F,l-,2+,3-,4 



<13> 4 



+4 1 



2(12)(23)(34)(41) 
<13) 2 



F1f C (si23,S234;W^,S23) 



+ \ F 4 ™ (s 2 3, S34; S234) + ^ F 4 ™ (s 34 , S14; S134) 



4iV r 



1 



(24) V 4 ( 24 ) 



F 4 ™ (S23,S345S234) 



-[42]Li(s 2 34,S23! 



N f N s 



(12)(41)(23)(34) 
(24)3 



1 

4(24) 



F1f(s 2 3,S34;S234) ~ [ 42 ] L 1 («234 , S 23 ) 
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(23)(41)[42] 2 / (14) (23) [42] f 



(24) V 
(12) (34) f , 
^(24T ( 23 ' 



-^3(S234,S23) 



+ 



(1 <-> 2), (3 <-> 4) 



(1«3) + (2h4) + (1«3),(2«4) 



(5.10) 



while the rational part is given by [48], 



R 4 (H,l-,2+,3-A + ) 



N f N s 
1- -d- + 



[24] 4 



+ 



iV c iVj 12[12][23][34][41] \s 24 s 12i 
(1 «- 2), (3 <- 4) 



f 523^34 _ 3 £23£34_\ 
S 24 / 



- 1 (y'Wfal 
+{(1<- 3)} + {(2<-4) 

5.4 The NMHV amplitude A$ (H, 1+ , 2~ , 3" , 4 



;i<-3),(2<-4) 



(5.11) 



By combining the results for the NMHV amplitudes given in sections |3| and |4| and the 
rational $ amplitude of [46] according to eq. fl2.4| ), we obtain the Higgs NMHV-amplitude, 
which is symmetric under the exchange (2 4). The finite cut-constructible contribution 



is, 



F 4 (tf,l + , 2- 3- 4-) 



5 234 



4(l| P/f |2](l| P/f |4][23][34] 
mPH\l? 



(34) 3 m^ 



+ 



+ 



2si34 (2 \p H 1 3] [34] [41] 2s 134 (1 |pif 1 2] (3 \p H 1 2] (41) 
1 ( (3|^|1] 4 , (24) 4 mf, 



4 5 i 2 4 V<3|p«|2](3|p H |4][21][41] (12)(14)(2|p H |3](4| P/f |3] 

ml(K\2f(M)* 



7=7± (PH">P1+P2) 



7 ( 7 + m2)(^l)(^3)(^4)(12) 



F| m (m^, Sl2 , S3 4) 



+ 1 



N f \f(3\p H \lf 



4NJ Vsi24[24] 2 
4(24)(3|p H |l 12 



F4F ( s 12) S14! S124) 



+ 



S124[42] 

' 7 

2s 124 (24)(34) 2 [41] 2 



f 4(23)(4|p H |l] 2 f 

M (S124, S12 J H r^TTvi ^1 ( s 123, «12 J 



si 23 [32] 

JV> JV.W [12][41](3|pg|2](3|p H |4] lm 
1 " iVe + iV^J V 2 S124 [24] 4 F4F (S12 ' Sl4 ' Sl24 ^ 



3 [42] 



L3 (S124,S1 2 ) 



, (34) [41] (3si 24 (34)[41] + (24) (3^ 1 1] [42]) f 

H iTFTTTTo ^2 («124, «12j 



+ 



3 [42] 2 

2s 124 (34) 2 [41] 2 (24)(3|ph|1] 2 



<24>[42]i 



3si 24 [42] 



^1 (S124,S12) 
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, (3|p H |l](45i24(34)[41] + {3\p H \l}(2s u + s 24 )) f , . 
+ Sl24 (24)[42]3 Lo (S124 ' Sl2) 

2 gl23 (23)(34) 2 [31]% (23)(34)[31](4|p H |l] f 
3[32] L3 (Sl23 ' Sl2) + 3[32] 2 (Sl23 ' Sl2) 

I (23)(4|p g |l] 2 f 
+ 3,s 123 [32] M-ia^ia) 



+ j(2~4)j. (5.12) 

For convenience we have introduced the following combinations of the finite pieces of one- 
mass (F^p) and two-mass hard (F^p h ) box functions (see Appendix |B|) , 

WW = F^(s 23 , S34; S234) + F^ h (s 4 i, s 234 ; "in, s 23 ) + Flf h (s 12 , s 234 ; S34, m 2 H ) 
= F^(si 4 , s 34 ; S134) + Flp lh (si 2 , si 34 ; m 2 H , s 34 ) + F^ h (s 23 , si 34 ; s u ,m 2 H ) 
W iz) = F^p 1 (si 2 ,si 4 ;si 24 ) + F^ h (s 23 ,si 24 ;m^,si 4 ) + F^ h (s 34 ,si24;si2,m^). 

In addition, to simplify the coefficients of the three-mass triangle F| m (K f, K^K 2 ) with 
three off-shell legs Kf, K%, K% 7^ 0, we use the notation of eq. ( |3.21| ). The rational 
part of the Higgs NMHV amplitude is given by eq. Q4.3| ) (which incorporates the rational 
A£\<tf,l + ,2r,3-,4r) amplitude derived in [46]), 

„,„ 1+rrn . !(, Nj, N 3 \l{ {23)(te)(4\p H \l][31] (3\p H \l} 2 
4 ' ' ' ' j " U N c + N c j2\ 3 Sl23 (12)[21][32] s 124 [42]2 

(24)(34)(3| P ff|l][41] _ [12] 2 (23) 2 _ (24)( S23 s 24 + S23 s 34 + s 24S34 ) 

3si 24 si 2 [42] Sl4 [42] 2 3(12) (14) [23] [34] [42] 

[2\vh\1](4\vh\1] 2[12l(23>[311 2 ' 



<2jggjjj(4jgHil] 2[12](23)[31] 2 ^ f 1 
' 3 S234 [23][34] " 3[23] 2 [41][34]J j + | (2 ~ 4) |' ^ 



6. Numerical Evaluation 

In this section we provide numerical values for the helicity amplitudes given in the previous 
section at a particular phase space point. To this end, we redefine the finite part of the 
Higgs amplitude as: 



-M 2 



A { P (H, l Al , 2 A2 , 3 A3 , 4 A4 ) = c r ^ (0) (H, l Al , 2 Aa , 3 As , 4 A4 ) ( - \ ^ 

+M% ' s (Ai, A 2 , A 3 , A 4 ) + ^Mf J (A!, A 2 , A 3 , A 4 ) + ^m(> s (\ 1 ,\ 2 , A 3 , A 4 )^j 

We evaluate the amplitudes at the phase space point used by Ellis et al. [19], 

p^ H = (-1.00000000000, 0.00000000000, 0.00000000000, 0.00000000000), 

= (+0.30674037867, -0.17738694693, -0.01664472021, -0.24969277974), 
p% = (+0.34445032281, +0.14635282800, -0.10707762397, +0.29285022975), 



(6.1) 
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p» = (+0.22091667641, +0.08911915938, +0.19733901856, +0.04380941793), 

p% = (+0.12789262211, -0.05808504045, -0.07361667438, -0.08696686795). (6.2) 



The results are presented in table |] where we have chosen the renormalisation scale to be 



Helicity 


A m 


Mi'* 








-116.526220-18.6817751 


-9.540396-0.0010101 


-0.176850+0.0010101 


0.176850-0.0010101 


H 


10.308088-0.824204 I 


-10.809925+0.0566461 


-0.388288+0.1983691 


0.296783-0.1551321 


h+ 


20.511457-0.888525 I 


-10.991033+0.3200091 


0.268501-0.0684141 


0.066595-0.015451 1 


- + -+ 


4.683784+4.242678 I 


-10.332320+0.1492161 


0.028668-0.066437 1 


0.166800+0.0388441 



Table 1: Numerical values for the finite parts of the Higgs + 4 gluon helicity amplitudes at the 
phase space point given in eq. (6.2). 



7. Conclusions 

We have calculated the last (analytically) unknown building block of the Higgs plus four 
gluon amplitude, the 0-NMHV amplitude 1 + , 2~, 3~, 4~). We chose to split the 

calculation into two parts, one being cut-constructible (to which we applied the techniques 
of four-dimensional unitarity) and a rational part, which is insensitive to four-dimensional 
cuts. 

By employing a four-dimensional unitarity-based strategy, we reconstructed the cut- 
constructible piece as a combination of bubble-, triangle- and box-functions. The coef- 
ficients of the box-integrals were obtained using quadruple cuts [28]. We used Forde's 
Laurent expansion method [29] to derive the coefficients of the one- and two-mass trian- 
gles. The only role of these functions in the amplitude is to ensure the correct infrared e 
pole structure, and as such they do not appear explicitly in our result. The coefficients of 
the three-mass triangle, and the bubble-functions were calculated using two independent 
techniques; the Laurent expansion method, [29] and the triple-cut spinor-integration tech- 
nique [30]. Finally the double cuts were computed analytically via Stokes' theorem [52] 
and checked numerically with the Laurent expansion technique [29]. 

To evaluate the remaining rational piece of the amplitude, we found it convenient to 
separate the rational part into two terms, one which is sensitive to the number of active light 
fermions Nj and a homogeneous part which is not. The Nf dependent contribution was 
efficiently computed using Feynman diagrams, while the homogeneous piece was amenable 
to a BCFW recursion relation approach [40, 41]. 

Previous calculations of 4> plus gluon amplitudes have considered c^-MHV, [47, 48], 
the 0-all-minus, [49] and the (/>-all-plus and (/>-nearly-all-plus rational amplitudes [46]. We 
have collected the results of these previous papers, so that the virtual one-loop amplitudes 

3 We have been informed by John Campbeif, that the entries for M^' 9 and M^' q in Table ^ are in 
agreement with results obtained using the seminumerical code described in Ref. [19]. 
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for the — > Hgggg process (in the large-mt limit) are available in one place. In partic- 
ular, we present compact formulae for the helicity configurations (H, 1 _ , 2~, 3~, 4~), 
A ( l ] (H, 1- , 2- , 3+ , 4+ ) , A { 1 ] (H, 1~ , 2+ , 3~ , 4+ ) and (H, 1+ , 2~ , 3" , 4" ) in section g know- 
ing that all other helicity configurations can be obtained by parity, or by permuting the 
gluon labels. 

The Higgs plus four gluon amplitudes have been calculated numerically in Ref. [19]. 
However, we envisage that the analytic formulae hereby presented can be used to provide 
a faster and more flexible analysis of Higgs phenomenology at the LHC. 
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A. Tree- level amplitudes 

Compact analytic forms for the tree-level amplitudes relevant for A^ \<fr, 1 + , 2~, 3 _ , 4~) 
have been computed using both MHV rules [18] and BCFW recursion [46] . 

The first main ingredients are the four simple, all- multiplicity results: 

4 )(<M + ,...,n+) = 0, (A.l) 
4°)(<M + ,...,r,...,n+) = 0, (A.2) 

4°)(6 1+, . . . ,r, . . . ,r , ■ ■ ■ , O = n »_ ( (l 4 + i) ' (A - 3) 
4Q)( ^ r '---' ra> IM ' (A - 4) 

In addition we need the non-trivial NMHV-type amplitude, 

4°)(<M + ,2-,3-,4-) = 

™J< 24 > 4 + (4M1] 3 (2|P |1] 3 



si24<12)(14)(2|p |3](4|p |3] si2 3 <4|p^|3][12][23] s lu (2\ P(f> \3] [14] [34] ' 
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The amplitudes for the <p coupling to a quark pair and a pair of massless scalars are 
given by, 



4°>(0, i s ,...,r,...,n s ) 
l- ...,r 
1+ ...,r 

with the non-trivial NMHV amplitudes given by, 



i • • • > n q 



i ■ ■ ■ ) n q 



(lz) 2 (m) 2 

(li) 3 (m) 

(-l) n (li)(m) 3 
nLi(^ + l)' 



4°)(^,l+,2-,3-,4r) 
(24) 3 m4 



<4|p |l] 



+ 



(2|^|1] 2 (2|^|4] 



5 i24<14)(2|p^|3](4|p |3] (4|p |3][12][23] s 13 4<2]^|3] [14] [34] ' 
4°)(</),l s ,2-,3-,4 s ) = 

si 23 (4b |l] (2|p |l](2|p |4] : 



(12) (24)^ 



+ 



Sl24 (14)(2|p |3](4|^|3] <4|p |3][12][23] Sl34 <2|^|3][14][34] 



(A.6) 
(A.7) 
(A.8) 



(A.9) 



(A.10) 



B. Scalar Basis Integrals 



P 2 
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7 



■p2mc 
r 4F 



Figure 4: Conventions for labelling the three scalar box integrals appearing in the one-loop Hgggg 
amplitudes. 

In this appendix we present the basis integral functions used to construct the various 
finite contributions to the Higgs helicity amplitudes. Figure^] sets our labelling conventions. 
The finite parts of the one-mass and two-mass (easy and hard) have the following forms, 



Fif (s,t;P z ) = -2(Li 2 ( 1- — ) +Li 2 ( I 



P v 



t 2 \t 6 



i2mh i 
4F I 



s,t;P 2 ,Q 2 ) = -2(Li 2 ( 1- — ) +Li 2 ( 1 - ^ ) + ^ ln z ( ~ 



P 



t 



l 



i ln (^) ln (^ 



(B.l) 



(B.2) 
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p2\ / p2\ / n2 

2me/ .. D 2 /~>2\ nl t : I i ^ \ , T : I i ^ \ r t ; / 1 



F^C*,*;^,^) = -2( Li 2 ( 1 — 1 +Li 2 ll- T ) +Li 2 ( 1 

+Li 2 [L-<£)-uJl- ^f-) + I In 2 (1) ) . (B.:!) 



t J \ st ) 2 \t 

The finite three-mass triangle is given by [58, 59], 

WM . ^± ( Li2 (- (i±f )) - U 2 (- , ( B„ 

where, 

3 

A = ^-M fc 4 + 2M|M| +1 (B.5) 

k=l 

Ml- Ml, x - ML. 2 / s 

4 = — = ^±4 (B.6) 

v A 



Alternative representations for these integrals can be found in references [60-62]. 
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